Mode-coupling theory and the 
fluctuation-dissipation theorem for nonlinear 
Langevin equations with multiplicative noise 

Kunimasa Miyazaki and David R. Reichman 

Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 
10027, U.S.A 

E-mail: km22399columbia.edu 

Abstract. In this letter, we develop a mode-coupling theory for a class of 
nonlinear Langevin equations with multiplicative noise using a field theoretic 
formalism. These equations are simplified models of realistic colloidal suspensions. 
We prove that the derived equations are consistent with the fluctuation-dissipation 
theorem. We also discuss the generalization of the result given here to real 
fluids, and the possible description of supercooled fluids in the aging regime. We 
demonstrate that the standard idealized mode-coupling theory is not consistent 
with the FDT in a strict field theoretic sense. 
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1. Introduction 

Mode-coupling theory (MCT) has been a useful first principles approach for studying 
the dynamics of supercooled liquids (see, for example, Ref.pQ). MCT was originally 
derived using projection operator methods together with several uncontrolled 
approximations. The theory has been successful in providing a quantitative account 
of many experimental and numerical observations, such as nonergodic parameters and 
structural relaxation exponents^- The approximate nature of the conventional MCT, 
however, has restricted its validity to description of only (i) two-point correlation 
functions^, (ii) for systems at equilibrium, and (iii) at relatively high temperatures 
below which the theory predicts a spurious glass transition. 

Experimental and numerical studies, on the other hand, have provided us with a 
rich dynamics, none of which conventional MCT can explain. For example, simulations 
and experiments of supercooled fluids have revealed the existence of the correlated 
local dynamical heterogeneities(see, for example, Ref.[5]), and dynamical scaling 6 . 
Another example is supercooled liquids brought out of equilibrium by quenching the 
system to low temperature or by adding shear flow. Here many experiments and 

| See, however, G. Biroli and J. -P. Bouchaud 3 • Note that, strictly speaking, the mode-coupling 
approach taken here is of the field-theoretic, diagrammatic variety. This is what allows the authors 
to compute multi-point correlation functions, something that is difficult to do in a useful manner 
from the projection operator approach. The field theory used in the work of Biroli and Bouchaud is 
that of Das and Mazenko|I] which assumes the FDT from the start. The work presented here is a 
first step beyond such a simplified treatment. 
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simulations have shown violation of the fluctuation-dissipation theorem (FDT) and 
the existence of the effective temperatures |S1 EJ- A microscopic theory which 
goes beyond MCT is desirable to account these phenomena. For such purposes, 
a systematic field theoretic approach is a good candidate. Such a field theoretical 
perturbation scheme for dynamical processes has been developed by Martin, Siggia, 
and Rose (MSR)JOJ- This approach is advantageous over the projection operator 
technique in that (i) it is conceptually straightforward to extend MCT-type equation 
to incorporate higher order moments of fluctuations, (ii) it is useful for calculating 
multipoint correlation functions which are essential observables to monitor the 
dynamical heterogeneities in the supercooled fluids, and (iii) it enables one to treat 
the correlation function and the response function (which are related by the FDT if 
equilibrium holds) on an equal footing and, therefore, it is a powerful tool for the 
treatment of noncquilibrium systems. 

Conventional MCT used for glassy or disordered systems is believed to be 
equivalent with renormalized perturbation theory [TT] without vertex corrections 
within the standard loop expansion of the MSR formalism|10|. This is true for a 
certain class of disordered systems such as the p-spin spin glass models [T2J - For 
supercooled fluids, however, this field theory has never been systematically used to 
derive the MCT equation even at equilibrium. The primary obstacle is that, although 
the original equations of motion satisfy the FDT at equilibrium, this does not imply 
that an arbitrary perturbation scheme also preserves the constraints of the FDT at 
each order of the expansion. Derivations of MCT for supercooled fluids from the field 
theoretic point of view has been discussed by several authors |%1 IT31 but either the 
FDT has been assumed (rather than consistently derived) ^ E] or a certain model 
has to be introduced for the derived equation to guarantee the FDTJ3]- Difficulties 
in deriving the MCT equations and extending them to higher order by a systematic 
loop expansion for supercooled fluids are due to certain properties of the nonlinearities 
in the microscopic Langevin equation. These difficulties do not exist in the schematic 
p-spin models. To illustrate the difficulties, let us consider a Langevin equation which 
describes dynamics of the density field, p(r, t), of the dense colloidal suspension, as an 
example 15 ; 

dp(r,t) 
dt 

where D is the diffusion coefficient which is assumed to be a constant, Sp(r, t) = 
p(r,t) — po with po = (p(r,t)) is a density fluctuation, f p (r,t) is a random noise, and 
c(r) is the direct correlation function. The second term in the brackets accounts for the 
interaction between the particles and is the source of the nonlinearity of the Langevin 
equation§. Eq.lQJ is known to exhibit the glassy properties at high densities. It has 
been also shown that, in equilibrium, this equation can be reduced to a standard 
MCT equation in the overdamped limit, using projection operators^! EH ED]- I n 
order to see the difficulty in applying the field-theoretic MSR, let us rewrite eq. QJ in 
a following form; 

+ «'.*>, (2) 

where 

Lp(v)p(r') = fciT'-DV • Vp(r, t)S(r - r') (3) 

§ Eq. is already coarse-grained in a sense that the bare interaction potential is replaced with the 
effective potential — c(r)/k-gT. See Refs. irp1ll7l 



DV-{vp(r,t)-p(r,i)V ^dr' c(r - r')Sp(r',t)\ + f p (r,t), 



(1) 
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is the Onsager coefficient and S is the entropy of the whole system which is given as 
a functional of the density by 

S = k B ^~ Jdr p(r)[bx{p[r)/p }-l] + ^Jdr J dr' c(|r - r'\)Sp(r)5p(r')j . (4) 

The random field f p (r,t) satisfies 

(f p (r, t)f p (r', t')) = 2k B L p{r)p{rl) S(t - t'). (5) 

One sees that there are two types of nonlinearities entangled in eq.J2J); One is due 
to the thermodynamic force SS/Sp(r), which is nonlinear in 8 p. Another is the 
density dependent Onsager coefficient which, through eq.JSJ, makes the random force 
a nonlinear function of the density, producing multiplicative noise. These properties 
are quite general for Langevin equations for realistic fluids. As elucidated in the next 
section, these two nonlinearities are the origin of difficulties which hamper the field 
theoretical approach for realistic fluids. 

In this letter, we develop the tools to treat some of the difficult nonlinearities 
discussed above, including multiplicative noise. Our goal is to derive MCT-type 
equations which satisfies the FDT at the lowest level of the loop expansion if the 
system is at equilibrium. This goal is an important prerequisite condition when we set 
out to explore nonequilibrium systems. Our theory will serve as a first step, however 
incremental, to prepare for the development of a field theory for supercooled fluids 
which goes beyond the conventional MCT and is capable of exploring, for example, 
nonequilibrium systems, the effect of higher order loops, and multipoint correlation 
functions. In particular, we will show that: 

(a) Care must be exercised in field theoretical derivations of the response of a system 
described by eq.JQ) to an external field since the response function is not trivially 
connected to the propagator in general. This fact has been ignored in past works pi], 

(b) The multiplicative noise term is essential for a proper treatment of the memory 
term. This point has been overlooked in past works [Tl]. as the proper field theoretic 
treatment of this term is subtle. 

(c) A standard one-loop treatment can lead to the usual MCT of Gotze and 
coworkers^ E] at the expense of satisfying the FDT. This will lead us to consider a 
slightly simpler model for which the associated self-consistent one-loop theory presents 
the FDT. This model will allow us to make connection with the work of Schmitz et 

a i.\n\. 

In the next section, we give more detailed accounts of the background and motivation 
of the present work. Section 3 is devoted to derivation of MCT for multiplicative noise 
using the MSR method. Consistency with the FDT is discussed in Section 4 and in 
Section 5. Perspectives for developing MCT for the nonequilibrium case are discussed 
in Section 6. 

2. Background 

In order to make the argument general, let us consider a classical stochastic dynamical 
process of a field variable Xi(t), where i is an index that denotes the type of field (such 
as the density) and coordinates (such as positions or wave vectors) which may be 
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either discrete or continuous. If the system is macroscopically at equilibrium, x,(£) 
obeys a nonlinear Langevin equation of a general form: 

dS T dS , „ dS 

Xi = M ia -K 1- L ia j: h Ji = A iQ h Ji, (6) 

where a sum over the repeated Greek indices is assumed. S is the entropy of the 
entire system. K^j = Mjj + is a kinetic coefficient which generally depends on x. 
M ia dS/dx a represents the reversible term where Af,j = —Mji is an antisymmetric 
matrix. L ia dS/dx a represents the irreversible term, where Ly is the (x-dependent) 
Onsager coefficient. /i(x, t) is a random noise which satisfies 

(/ 2 (x,t)/,(x',t')) x(t)=x = 2fcB^(x)5(i-0, (7) 

where (■ • -) x ( ( ) =x denotes the conditional average in which the ensemble average 
is taken with a fixed value of x(i) = x at time t. The fact that the Onsager 
coefficient is a function of x means that the random noise is also the function of 
x, i.e., it is multiplicative [22] . Eq.© is the general expression for dynamical systems 
whose stationary distribution function in the absence of nonequilibrium constraints is 
given by the equilibrium ensemble. Examples include the (fluctuating) Navier-Stokes 
equation |2H1 E] and the nonlinear diffusion equation (eq.QJ). The response function 
Xij(t) is defined as the response to a time-dependent external force F(t) by 

(A Xi (t)) F = f dt' xUt - t')F a (t'), (8) 

J — oo 

where (Axi(i))f is the deviation of (xi(t)} from its equilibrium value due to the 
external force. The FDT asserts that Xijif) is related to the correlation function in 
the absence of F(t), Cij(t) = (xi(t)xj(0)), by 

*«<"> = -^fT fat ^- < 9 > 

The FDT is proved easily using linear response theory [2SJ - The FDT is one of the 
strongest and the most robust statements of equilibrium statistical physics and it holds 
for any dynamical processes, classical or quantum, linear or nonlinear, as long as the 
system is stationary, satisfies the condition of detailed balance, and the perturbation 
is small enough. 

The MSR formalism allows us for the use of a systematic loop expansion for the 
solution of the nonlinear Langevin equation in terms of the moments. MCT is regarded 
as the lowest order self-consistent approximation with no vertex correction in the MSR 
formalism. Deker et al. [5T] have proven that the FDT holds at each order of the loop 
expansion for three special classes of dynamical processes: "Class A" where, in eq.©, 
My(x) = 0, Lij is a constant (independent of x) and thus the noise is additive. The 
nonlinearity of the Langevin equation originates from the entropy S. "Class B" where 
L,^ is constant and the entropy is a quadratic function of x. The reversible matrix 
Mij (x) depends on x which is the origin of the nonlinearity. The equations discussed 
by Kawasaki to describe dynamical critical phenomena[2H| belong to this class. "Class 
C" involves Hamiltonian systems which do not have an irreversible part. 

The problems is that, as discussed in the Introduction, even eq.([Q), which 
describes the dynamics of supercooled fluids, does not belong to any of the classes 
listed above. The nonlinear term in eq.([Q) originates from the combination of the non- 
quadratic term of the entropy and the variable dependence of the Onsager coefficient. 
Extensions of the MSR formalism to more general cases have been discussed in 
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Refs. |14l 1271 1281 I29| but derivations of MCT equations have not been given so far. 
Since Deker's classifications do not cover these dynamical processes, it is convenient 
to re-categorize the nonlinear stochastic processes in a slightly more general way than 
Classes A-C of Deker: 

(Class I): The nonlinearity comes solely from the entropy. Kij is independent of x. 
The noise is additive. The mean-field model of p-spin spin glasses also belongs to this 
tvpe|12j. Again, it is trivial to show that the FDT holds at each order of the loop 
expansion. 

(Class II): The entropy is a quadratic function of x but Kij(x) is dependent on x. 
Lij(x) can be also a function of x and, therefore, the noise can be multiplicative. 

(Class III): The entropy is an arbitrary function of x and Kij(x) is dependent on x. 
Real fluids including the one described by eq.Q belong to this type. 

In this letter, we shall focus here on Class II. Here we restrict ourselves to the simplest 
situation where K^ (x) is a linear function of Sx: 

K iJ (x) = Kg ) +Kg! a 5x a , (10) 

because this is the most important case in the context of the glass transition. We shall 
show, for this class of Langevin equation, that MCT is consistent with the FDT even 
with the presence of multiplicative noise. We believe that this conclusion is valid for 
arbitrary function of Kij (x) of Class II. We will end this work with some comments 
on the more interesting Class III case which is problematic from the standpoint of the 
FDT within the MCT approximation. 



3. MSR formalism for processes with multiplicative noise 

In this section, we shall develop the MSR method for the Class II case. In Class II, 
the entropy in eq.l© is given by a quadratic form 

S = S + -tonpSXaSxp, (11) 

where Sx = x — (x) and fi^- = —knC^{t = 0) is the inverse of the equal time 
correlation function. Such an approximation is, perhaps, not as crude as it may 
superficially appear. Indeed, it has been shown via direct simulation that both 
simple[3U] and complex liquids]!^ have Gaussian density fluctuations over a wide 
range of length scales. Ea.ljlip. within the canonical ensemble, is a precise statement of 
this approximation. The Fokker-Planck equation for the probability density function 
P(x, t) equivalent with eq.JBJl is written as 



8P(x,t) d r , x dS , T , . d \ „. . , 



(12) 



For the detailed balance condition to be satisfied, Mj,'(x) must satisfy the following 
condition (potential condition|25j): dMi a (x)/dx a — 0. We shall also assume the 
similar incompressible condition for L !3 (x): dLi a (x) /dx a = 0. These two conditions 
are satisfied for most hydrodynamic equations including eq.QJ (see eq.Q). The latter 
condition is especially useful because it enables us to avoid distinguishing between 
the Ito and Stratonovich interpretations which are associated with the multiplicative 
noise(22]- Using eas. ljlQ(l and Ijlljl . eq.® can be rewritten as 

Xi = fii a X a + -ViapXaX/j + fi. (13) 
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Hereafter we shall omit the lt S" in front of x. In eg. l|13|) . My — K^Cl a j is a bare 
transport coefficient and V«n, = Vikj is the symmetrized vertex defined by 

v m = M^n ak + Mj£ k n aj + iflp ak + Lfl k n aj = M ijk + c ijk , (u) 

where Aiijk and C^k are reversible and irreversible contributions of Vyfc, respectively. 
Mijk satisfies the cyclic condition given by 

^iaM a jk + QjaMaki + ^k a M mJ = 0. (15) 

This is proved using the condition that the reversible part does not contribute to the 
entropy production. 

Following the standard MSR procedure|l()|. we shall introduce a spinor z = (x, x), 
where Xi = —d/dxi. It is convenient to define the generating function for z by 



W[$] = In ( exp 



dt £ ■ z 



(16) 

where £ = (r/,77) is the auxiliary field conjugate to z which will be eventually set to 
zero. "exp + " implies the time ordering which aligns the quantities with larger t on 
the left. We define the cumulant function (1, • • • , n) = ((z(l) ■ ■ ■ z(n))) by 

5 n W\£] 

o^-^-mA- (17) 

The index number 1 = (i,t, ±), etc... represents the index for field variables i, time t, 
and the index of the spinor defined by z(+) = x and z(— ) = x, respectively. Let us 
construct the equation of motion for the first cumulant, ((z(l))). Substituting ca. (|13fl 
into ea. (|17fl for n = 1, we obtain the Schwinger equation: 

= (([T, *(1)])> + ia(l, l)e(l), (18) 

where a sum over the repeated underlined indices is assumed. [■ • •] is the commutator, 
T is the Fokker-Planck operator defined by ea. l|12|) represented in terms of (x, x), and 

ia(l,2)= ( ° - 1 )8(h-t 2 ). (19) 





The explicit expression of ea. l|18|l is given by 
d((xi}) 1 

dt =f)i + Via ((X a )) + -V iaf } {((x a Xp)) + ((x a )) ((xp))} 



2k T 



4°J «**» + L^p {((x a Xp)) + «£ a » ((Xp))} 



d 

dt 



(20) 



-T)i- Via ((x a )) - V a pi{((x a Xp)) + ((x a )) ((xp))} + V aai 

-kBL%{((x a xp)) + ((x a ))((xp))}, 

where '/Xjj is the transverse of fiij. The last terms in these equations are due to the 
multiplicative noise. Ea. H2Q(l is written in short as 

So-'a.U «*(!)» = S(l) + C(l) + i 7 s(l, 1,2) {((z(l)z(2))) + «*©»} , (21) 

where C(l) = (V aa i,0) is a constant which does not contribute to the following 
arguments and Qo(l, 2) is the bare propagator whose inverse is written as 

g - 1 (l,2)=za(l,2)^+ 72 (l,2) (22) 
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with the symmetric matrix 72 (1, 2) defined by 

72(1,2) =(_° M _ 2 -^ (0 ))^-* 2 ). (23) 

The non-zero components of 73(1, 2, 3) are 

[lz{h,h,-;i2,h,+;hM,+) = Vi li2 i 3 5(ti -t 2 )S(t 1 -t 3 ) 

{ (i) ( 24 ) 

[73(11, t\, -;i2,*2, -;«3,*3,+) = 2fc B L- 1 - 2 , l 3 (5 ^ 1 -*a)<K*i -*s) 

and its permutation of the indices (1,2,3). Note that 73(1,2,3) is a fully symmetric 
tensor. The second moment £(1,2) = £?£ = o(l,2) = (z(\)z(2)) is given by taking the 
derivative of ea. 1)21(1 with respect to ((z(2))) using ea. 1)17(1 and then turning off £ = 0: 

£- 1 (l,2) = £ - 1 (l,2)-S(l,2), (25) 

where S(l,2) is the self energy. By neglecting the vertex correction, we obtain the 
MCT expression for the self-energy: 

£(1,2) ~ ^73(1,1,2)£(1,3)£(2 ) 4)7 3 (2,3 ! 4). (26) 
We can write the components of these matrices as 

=(g t and E(l,2)=(2 g ) , (27) 

where "f" represents Hermitian conjugate defined by — £') = — t). From 

the structure of ea. f2"5|) . it is straightforward to show that </(—,—) = £(+,+) = 0. 
Cy(i — t') = (xi(t)xj(t')) is the correlation function and Gij(t — t') = (xi(t)xj(t')) is 
the propagator which describes the response of the system to the random noise. For 
t > 0, the equations for C(t) and G(t) can be written explicitly using eas. ((22(1 - 1(25(1 as 

dC (t) /"* 

— ^ =Mi«Caj(i)+ / dti B<a(i-ti)C^(ti)+ / dti Dfa(t-*i)G^(ii) (28) 



— OO 



(30) 



d/ fJ,iaG aj (t)+ J dh Eiait-t^Gajftl) (29) 

with the self energies given by 

Eij(t) =Vi a pG a \(t)Cfj ti (t)V\^j + k^ViapGaX^Gfj^t)^^- 

Da® =\v iaP C aX {t)C Pll {t)V 3 x li + 2k B V tal3 G aX (t)C l3t ,(t)L { j L l^ 

where use has been made of the causality condition: G(t) = for t < 0. The terms 
containing in ea. ((30(l originate from the multiplicative noise. Similar terms were 
derived by Kawasaki et al.\TQ but they were disregarded and their importance was 
not addressed. Note that the propagator G(t) represents the response to the noise but 
it is not the response to the external force defined by eq. 0). The response function is 
obtained by evaluating the linear response of the average (xi(t)) to the external force 
F(t). The term associated with the external force is introduced naturally by replacing 
the entropy with the one associated with the work done by the force as 

S F = S+^. (31) 



MCT and the fluctuation-dissipation theorem 



8 



Inserting this expression into the entropy term in eq. and taking the leading order 
of the formal solution, it is straightforward to derive the expression for the response 
function. The result is 

XV (t - = ^{xi{t)x a {t')K aj {x.{t'))) 

= -( Xl (t)£ a (t'))K^ + ¥ (x t (t)£ a (t') X0 (t'))K^. 

The three point correlation function (xi(t)x a (t')x/3(t')) in this expression is calculated 
in the same spirit as derivation of ea. (|26[l . Up to the one- loop level, neglecting 
the vertex correction, it is written as (z{\)z{2)z{2>)) ~ (5(1, 1)£(2, 2J£(3, 3)73(1, 2, 3). 
Substituting ea. (|24[) to this, one obtains the MCT expression of the response function: 

Xij(t) = ^G ia (t)K^ + i J*dh G ia (t - hp^G^C^K^. (33) 

In previous works, the propag ator T- 1 G(t) • has been called the response 

function[2J. But as discussed above, it is not identical to the full response function 
in general. They become identical only if the kinetic coefficient is a constant, i.e., 
for Class I. 



4. FDT and MCT for Class II 

In this section, we prove that the correlation and response functions given by eas. H28|) . 
((23, and (J23) satisfy the FDT, eq.Q. Deker et al. have shown that, for the MCT 
equation of Class B, there is a simple relation between the correlation function and 
the propagator if Ly is a constant or — 0|21|: 

G(t) = e(t)C(t) ■ C- x (0), (34) 

where 9(t) is the Heaviside function. We prove that this is also true when ^ 
or the noise is multiplicative as follows: Taking the time derivative of both side of 
ea. (|34|) and substituting the equation for the correlation function, ea. (|28fl . we have 

— = 1 + (j, ■ G + 9(t) {E <g> C + D ® G f } • C-^O), (35) 

where A <g> B = dti Ai a (t — ti)B a j(ti). The terms containing Aiijk in 

E (g> C ■ C _1 (0) can be rearranged using eq.(P3} and the cyclic condition, ea. lfT5|l . 
For example, if t\ < 0, 

V« tt ^ aA (r)^(r)M v a j (f 1 ) = —ViafiCaxWCfrWMuXvGtljfa), (36) 

where r = t — t\. Here we have used the fact that CV,-(0) = — fcefi"- 1 . On the other 
hand, the terms containing djk in E <g> C • C _1 (0) are rearranged as 



^Gax^Cp^^Cx^k + kn 

= {-2k B G aX (T)C^(T)L[^ - l -C aX {T)C^{r)C kXl ^ G^h). 



(37) 



Eq.||2ni) combined with eq.(|37|) cancels with D (g • C^ 1 (0) of eq. (JSHJ - Likewise, 
for t\ > 0, E <S> C • C _1 (0) can be rewritten as E <E> G^. Therefore, eq. l|35|) becomes 
equivalent to the equation for G, ea. (|29|l . This is the end of the proof. 
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Now let us prove the FDT. By taking the derivative of ea. (|34[l with respect to 
time and using equation for G(t), ea. l|29[l . we have 



dC(t) 



-k B G(t) • K (0) + f dtx G(t - h) • E(ii) • C(0). 
Jo 



(38) 



dt 

In this expression, E(£) • C(0) can be rewritten using ea. H34|) and antisymmetric 
property of My \ as 

{E(t) • C(0)}y = - k B V ia pG aX (t)C^(t)K^. (39) 

Therefore, the right hand side of ea. (|38|l becomes identical to —ksTxijit) given by 
eq. <|33|) . Thus we arrive at eq.@ and the FDT is proved. 

Finally let us derive the closed equation for Cy (t). D ® again cancels with 
E ® C for t\ < 0. For ti > 0, Eij(t) is rewritten as 

En® = -\v ia pC aX (t)C f3 ^t) (V - 2C) vXll (7-/(0) (40) 
and we arrive at the MCT equation for Cy (t) ; 



^ a C aj (t) + / dti S m (t - *i)C aj -(ti). (41) 



o 

It is important to realize that —2C in the vertex in eg. \4U\j is due to the multiplicative 
noise and the presence of it is essential. For example, for the pure dissipative case 
(Mjj(x) = 0), neglect of the multiplicative noise leads to the wrong sign in front of 
the integral term (and thus leads to pathological behavior). This term is neglected in 
Ref.[TI|. 

5. FDT and MCT for Class III 

In this section, we shall consider the category of problems we call Class III and 
elucidate generic reason why the MCT approximation for the Class III dynamics is 
inconsistent with the FDT. For complete discussion of the technical aspects involved 
would require a longer discussion than we provide here. Our main point is simply to 
sketch the difficulties that arise in attempting to formulate a simple MCT (namely a 
self-consistent one-loop theory for both the propagation and the response function) 
that satisfies the FDT. The important conclusion is that the standard idealized MCT of 
Gotze and coworkers lj[2] can not be consistently derived via field-theoretic techniques, 
at least via the usual one-loop approximations applied to eq.Q. This is discussed both 
in this section and in Section 6. 

For problems of Class III, the entropy is not a quadratic function but has higher 
order terms; 

S = S + -Sl a p5x a 5xp + —AapySxaSxpSxy + ■ ■ ■ . (42) 

We again assume that the kinetic coefficient ify (x) is a linear function of x. Up to 
the quadratic order in x, the nonlinear Langevin equation for x is given by ea. (|13f) 
but the vertex, ea. H14|) . is now modified as 

= v"> + v<'<> <43) 
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where vS = K^-^-ajk is the vertex that originates from the nonlinearity of the 

entropy whereas V^P , which is the same as ea. (|14f) . originates from the x-dependence 
of the kinetic coefficient Kij(x). The one-loop equations for Cij(t) and Gij(t), eqs.(J5SJ 
and l(2"5|) , remain the same except that the vertex is given by eq. lj4*5|) instead of eq. lfH|) . 
First, we shall show that there is no simple relation such as ea. (|34f> which relates the 
correlation function to the propagator. The starting point is the formal solution of 
eq.J2Hl for the correlation function; 

C = G(g> (fc B K (0) + /c B K (0)t +D) ® G*. (44) 

This is derived by eliminating E(i) from the formal solution of C(t) by substituting 
the formal solution for G(i). Substituting the equation for G, eq.J55J), eq.(@3J can be 
rewritten as 

C = -kvG-n- 1 + G®f ® Gt (45) 

for t > 0. In this expression, f(t) = -fc B E(<) ■ - fe^T 1 ■ E(i) + D(t). Following 
the similar steps as eas. (|37tl3"9"|) . f(t) can be reduced as 

fij(t) = - V ia pG aX {t)C^{t)V { ^-j + l -V taP C aX {t)C^{t)vfl 
+ (Higher order loops), 

where only the first two terms (one loop) are explicitly shown. The other terms 
consist of the higher order loops. These loops appear to contain at least one multiple 
fV (/) V (//) . This means that these higher order loops do not appear in either Class 
I or Class II problems. If one takes the time derivative of ea. l)45|l . one obtains, after 
straightforward but tedious calculations, the following expression; 

= -k B Tx(t) + (Higher order loops). (47) 

dt 

Again, the higher order loops always contain at least one multiple of V (/) V("). One 
of the lowest order terms is an integral such as 

dt' ViapGaxQ - t')C fj ,{t - i')vSvS,^, A ,(t')^^(i')V iA ,^. (48) 

Note that this term is an irreducible loop in the field theoretic language, which means 
that this can not be represented by any simpler renormalized diagram. Since the 
original Langevin equation, eci. (|13fl . itself does satisfy the FDT, the failure of the 
FDT in ca. (|47|l is attributed to the inconsistencies of the loop expansion with the 
FDT. In other words, a naive loop expansion using the bare fields z = (x, x) for Class 
III problems does not preserve the FDT at the each level of expansion; The higher 
order diagrams shown in ea. (|47|l are cancelled only by taking the next higher order 
loops in the MSR loop expansion in Section 3. 

The failure to derive the FDT at the one-loop level for Class III leads to the 
failure of deriving a MCT type equation such as ea. 141(1 . By substituting eas. ((45|l and 
H47J1 into the equation for dj(t), ea. 1)28(1 . we obtain 

dC«j(*) 
dt 

where Eij(t) is given by the same expression as eq.(@UJ> except that V is given by 
ea.ljl3"|l and C in the second vertex, (V — 2£) in eq.|@UJ| is replaced by 



Hi a C a j(t) + / dti E ia (t — ti)C a j(ti) + (Higher order loops), 
Jo 



(49) 



L^Kjk + L^Mak + L^Aj. (50) 
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The higher order loops in eq. I|49|l are again irreducible diagrams which do not appear if 
either or V^ 11 * 1 is absent. Eq. l|49|) shows that the standard MCT equation for the 
Class III is regarded as the uncontrollable approximation in the held theoretic sense 
in that it neglects terms which are present in the original set of equations, eas. l|28|l 
and H29fl . The consequences of this are discussed below. 



6. Discussion 



In order to see how the result shown in the previous sections are related to real fluids, 
let us consider the Langevin equation for a colloidal suspension given by eq.Q. The 
Class II equation is derived by approximating the entropy given by eq.Q with its 
Gaussian form. Neglecting the terms of higher order than quadratic, one has 



2 



is the Fourier transform of Sp(r), N is the total number of the particles, and 
S(k) = ./V -1 (l^pkl 2 ) is the static structure factor. As mentioned in Section 3, this 
approximation has been shown to hold over a wide range of length scales in real 
liquids 30] |^. It is also compatible with the functionals used to derive integral 
equations of fluid structure [32] • It is not expected to hold on very short length scales 
where density fluctuations are Poissonian and the ideal gas entropy is essential. On the 
other hand, density fluctuations on such length scales are not expected physically to be 
effective in providing glassy behavior. It is plausible that the approximation ca. (|51|l 
may be used (along with an appropriate large wavevector cutoff) in the treatment 
of realistic fluids. Indeed, the approximation, ea. (|51|) . is used in Ref.|13|. Using 
eqs.(|3]) and (|51|l . the MCT equation, ea. l|41|l . for the density correlation function 
F(k,t) = N^ 1 (6pk(t)6p^ k (0)) is written as 

= ~W) F{kl t] + f dh M ' {k ' 1 ~ tl)F{k ' tl] (52) 

with the memory kernel given by 



o 



/>-'//-' /' dq fkq k-pV 



M '^ = 2^Jwf\W) + W)S F{q > t)W) > (53) 

where p = k q. Note that the vertex k • q/S(q) + k • p/S(p) is precisely the one that 
appears in Ref.|13j. where the Gaussian approximation to the entropy is also made. 
This equation should be compared with the standard MCT equation which has been 
derived using the projection operator method with the decoupling approximation^] 
where || 

M(k,t) = ^^- /-^-(k-qc^+k-pc^)} 2 F(q,t)F(p,t). (54) 



We observe that there is a difference between eas. H53|l and Ij54(l : l/S(q) appears in the 
vertex function of ea. (|53|l . whereas the direct correlation function poc(q) = 1 — 1/ S(q) 

|| Note that eg. 1521 is different from the MCT equation used in the supercooled fluids; —Dk 2 F(k,ti) 
appears in the place of 8F(k. fi )/dti IT51 1201 The difference originates from the overdamped 
nature of the starting diffusion equation, eq.0. One obtains the dF(k,ti)/dt\ term if one 
incorporates the momentum density as well as the number density as stochastic variables. Technically, 
this is equivalent to using the irreducible projection operator introduced by Cichocki et aL llHl . 
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shows up in ea. H54|) . This difference traces back to the Gaussian approximation, 
eq. l(5T|) 34 j . The entropy of fluids is not Gaussian due primarily to ideal gas part, p In p, 
in eq.Q). As discussed in Introduction, the nonlinearities of realistic fluids come both 
from the entropy (in this case, plnp) and the kinetic coefficient and therefore realistic 
fluids are destined to belong to Class III over the entire range of wavevectors. Indeed, 
if the full expression for S, eq.Q, instead of approximated form of ea. (|51|) is used, one 
sees that the non-quadratic term of S gives a vertex of the form of — D(k • q + k • p) 
which, combined with D{k ■ q/S(q) + k • p/S(p)}, leads to poD{k ■ qc(q) + k • pc(p)}. 
Using the full expression for S means that the dynamics now belongs to Class III. As 
illustrated in Section 5, one obtains the memory kernel in the form of eq.JMJ but there 
are always extra terms which are a direct consequence of inconsistencies of the FDT 
with the one-loop approximation in the MSR formalism, at least if the loop expansion 
is made directly with physical density modes as field variables^. This argument is 
true for arbitrary orders of the loop expansion. This conclusion implies that there 
is no simple systematic way to derive ea. <|54[) from cq.Q using the standard field 
theoretic method. This inconsistency is not relevant as far as one is concerned only 
with the equilibrium state, because one may always adopt an approximation where 
one neglects higher order terms in ea. H49(l and "define" the response function via the 
FDT instead of solving the equation for G(t), ea. i|29J) . separately. But we can not do 
so if the system is in nonequilibrium state. It is desirable to develop such an expansion 
method that preserves the FDT relation at each level of the perturbative expansion. 

Even with the difficulties discussed in this letter, the MSR or field theoretic 
approach is still an attractive route to attack out-of-equilibrium supercooled fluids, in 
that it is systematic and one does not need to evaluate the nonequilibrium measure 
which is required in alternative approaches such as the projection operator technique. 
Models which belong to Class I have been already discussed extensively in the context 
of spin glasses and even for supercooled fluids [SSJ. However, it is difficult to construct 
realistic models of Class I which can incorporate the effect of changes of structure 
embodied in S(k) observed in simulations of aging0 and sheared systems [Ml 137] . 
which is argued to be essential for the violation of the FDT of real fluids|3H|. The 
Class II system derived here is a better candidate as a realistic "model fluid" . As 
discussed above, the Gaussian approximation for the entropy is known to be a good 
description for wide ranges of densities and length scales PHI EI] The equation 
for equal-time correlation functions such as S(k) should be constructed in the same 
manner as the MCT equation derived here. The solution could be plugged into the 
vertex functions of the set of the MCT equations, eas. l|28l) and (|29|l . Eventually, 
these three equations can be solved self-consistently. Performing such a calculation 
will require some consideration of a wavevector cutoff that would eliminate spurious 
divergence that arise from the approximate vertex functions (ea. f53f> ). Calculations in 
this direction are underway. 

In summary, in the present work, we have broadened the range of applicability 
of the MSR approach by extending the method to dynamical processes with the 
multiplicative noise. This is a necessary step in the treatment of the Brownian 

% Recent unpublished work by G. Biroli, A. Lefevre, and J. -P. Bouchaud shows that if a transformed 
set of modes is used, FDT can be recovered for the full class III problem, although, at the time 
of submission of this work, other mathematical difficulties arise that render the resulting equations 
unusable in the present form. In particular, the vertex that results is distinct from that of eg. 1541 . 
and thus the one-loop approximation does not yield the usual form of the standard idealized MCT 
of Gotze and coworkers. 
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dynamics of colloidal suspensions. This formalism still does not cover real fluids 
over all length scales and it is found that there is no direct compatibility between the 
MCT equations derived from MSR approach and from the the projection operator 
method. We proposed an approximate but feasible method to explore nonequilibrium 
supercooled fluids using the formalism discussed in this paper. The MSR method is not 
restricted to evaluation of the two point correlation functions nor to the lowest order 
loop expansion. Extension of the method to the multipoint correlation functions and 
to higher order loops will be essential for understanding growing length scales which 
are hidden in the supercooled fluids : 3]. The formulation presented in this paper will 
serve as first step towards such extensions. 
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